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Abstract The Warm-Hot Intergalactic Medium (WHIM) is thought to contribute 
about 40 — 50 % to the baryonic budget at the present evolution stage of the universe. 
The observed large scale structure is likely to be due to gravitational growth of density 
fluctuations in the post-inflation era. The evolving cosmic web is governed by non-linear 
gravitational growth of the initially weak density fluctuations in the dark energy dom- 
inated cosmology. Non-linear structure formation, accretion and merging processes, 
star forming and AGN activity produce gas shocks in the WHIM. Shock waves are 
converting a fraction of the gravitation power to thermal and non-thermal emission 
of baryonic/leptonic matter. They provide the most likely way to power the luminous 
matter in the WHIM. The plasma shocks in the WHIM are expected to be coUisionless. 
CoUisionless shocks produce a highly non-equilibrium state with anisotropic temper- 
atures and a large differences in ion and electron temperatures. We discuss the ion 
and electron heating by the coUisionless shocks and then review the plasma processes 
responsible for the Coulomb equilibration and coUisional ionisation equilibrium of oxy- 
gen ions in the WHIM. MHD-turbulence produced by the strong coUisionless shocks 
could provide a sizeable non-thermal contribution to the observed Doppler parameter 
of the UV line spectra of the WHIM. 
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1 Introduction 



Cosmological simulations of tiie large-scale structure (LSS) predict that about 40—50 % 
of baryons at epoch 2 < 2 could reside in the Warm-Hot Intergala ctic Medium (WHIM) 
with temperatures 10^ — 10^ K at moderate overdensities S < 100 ( Cen fc Qstrikeill999l : 



bave et al.ll200ll : lFang et al.ll2002h . The WHIM heating is due to shocks dri ven by grav- 
itationally accelerated flows in the LSS structure formation scenario (e.g. iKang et al.l 
120071 '). Numerical simulations predict the observational signatures of the web g 
function of redshift. The simulations account for feedback interactions between galax- 
ies and the intergalactic medium, and demonstrate that the X-ray and ultraviolet 
Ovi, Ovil and Ovill lines and the Hi Lyma n alpha line are good tracers of low- 
density cosmic web filamentary struct ures (e.g. iHellsten et al.lll998l : iTripp et al.ll200Gl : 
ICen et ahlboOlMFurlanetto et al.ll2004l ). Intervening metal absorption systems of highly 
ionised C, N, O, Ne in the soft X-ray spectra of bright Active Galactic Nuclei (AGN) 
were suggested to be tracer of the WHIM. The predicted distribution of ion col- 
umn densities in the WHIM absorbers is steep enough to provide only a f ew systems 
with A^OviI ^ ^^^^ cm^^ along an arbitrary chosen line of sight (e.g. iFang et al.l 
Therefore, the detection of the WHIM is particularly difficult and requires 
very sensiti ve UV and X-ray detectors, both for absorption and for emission pro- 
cesse s (e.g^ Lehner et al. 20071: Nicastro et al. 20051 : iKaastra et al.l I2OO6I : iTakei et al.l 
|2007| and lmchter et al.ll2008l : lOurret et al.ll2008l - Chapters 3 and 4, this volume). Fu- 



ture projects and namely Cosmic Origin Spectrograph (COS), the X-Ray Evolving 
Universe Spectrometer {XEUS), Constellation-X and the Diffuse fntergalactic Oxygen 
Surveyor (DfOS) will increase the signal-to-noise ratio in the spectra allowing to study 
weak systems with A'jjj < 10^^'^ cm~^ and A^Qyjj < 10^^ cm~^. Simulations of 
spectra of the broad Lya absorption lines and the highly ionised oxygen lines in the 
weak systems require thorough modelling of physical con dition in the plasma (see e.g. 
lMewdll990l : |Paerels fc Kahnll2003l : lKawahara et allbood ). We discuss in this paper the 



heating and equilibration processes in the shocked WHIM plasma affecting the spectral 
simulations. A discussion of c oUisionless shock p hysics relevant to cosmological shocks 
in the WHIM can be found in lBvkov et al.ll2008l - Chapter 7, this volume. To this end, 
in this paper we discuss first some specific features of coUisionless shock heating of ions 
of different charge states providing highly non-equilibrium initial states just behind the 
magnetic ramp region that relaxes to an equilibrium state through Coulomb collisions, 
and the relation of the processes to simulations of emission /absorption spectra of the 
WHIM and observational data analysis. 



2 WHIM heating and ion temperature evolution 

The plasma ion heating in the WHIM is most likely due to cosmological shocks. The 
Alfven Mach number of a shock propagating through an ionised gas of local overdensity 
S — p/ (p) at the epoch z in the standard ylCDM cosmology is determined by 

Xa = Wsh(47rpi)^/Vs ~ 20.6vs7S^/^{l + zf/^{n^h^/0.02)^/^BZl. (1) 

where -B-g is the magnetic field just before the shock, measured in nC and Vgj is the 
shock velocity in 10^ cms~^, and (p) is the average density in the Universe. 
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The sonic Mach number for a shock propagating in a plasma of standard cosmic 
abundance is 

-8.5 i;s7 [14(1 + /ci)]"^/^ (2) 

where T4 is the plasma ion temperature measured in 10"^ K (typical for a preshock 
photoionised plasma) and /ei = Tc/T-^. An important plasma parameter is 

p = Mil Ml ~ 6<5(i + zf^n^y? iQm)Bzl\Ti(\ + /ci)]. 

It is the ratio of the thermal and magnetic pressures. In hot X-ray clusters of galaxies 
the beta parameter is ~ 100 for ~ [iG magnetic fields in the clusters. The most uncer- 
tain parameter is the magnetic field value in the WHIM allowing for both /3 ~ 1 and 
/3 ^ 1 cases. 

In a supercritical coUisionless shock the conversion of kinetic energy of an initially 
cold fiow to the ion distribution with high kinetic temperature occurs in the thin ion 
viscous jump. The width of the ion viscous jump Zivi in a coUisionless shock propagating 
through a plasma with /3 ~ 1 is typically of the order of a ten to a hundred times of 
the ion inertial length defined as = c/tjpi ~ 2.3 x lO^n^'''^ cm. Here tjpi is the ion 
plasma frequency. The ion inertial length in the WHIM can be estimated as 

h ~ 5.1 X 10^°5"^/^(1 + 2)"^/^(J7b/iV0.02)~^/2 cm. 

The width of the shock transition region for magnetic field is also > lOZ; for a quasi- 
perpendicular shock, but it is often about ten times wider for quasi-parallel shocks. 

Properties of nonrelativistic shocks in a hot, low magnetised plasma with high 
/3 S> 1 are ye t poorly studied. M easurements from the ISEEl and ISEE2 spacecrafts 
were used bv lFarris et al.l 1 19921 ') to examine the terrestrial bow shock under high beta 



conditions. These measurements were compared with and found to be in agreement 
with the predicted values of the Rankine-Hugoniot relations using the simple adiabatic 
approximation and a ratio of specific heats, 7, of 5/3. Large magnetic field and den- 
sity fluctuations were observed, but average downstream plasma conditions well away 
from the shock were relatively steady, near the predicted Rankine-Hugoniot values. 
The magnetic disturbances persisted well downstream and a hot, dense ion beam was 
detected leaking from the downstream region of the shock. The observation proved the 
existence of coUisionless shocks in high beta plasma, but a detailed study of high beta 
shock structure is needed for cosmological plasmas. 

We discuss in the next section the ion heating in coUisionless shocks illustrating the 
most important features of the process with the results of a hybrid simu lation of the 
oxyg en ions heating in a quasi-perpendicular shock considered earlier bv lBvkov et al.l 



boOSi - Chapter 7, this volume. 



2.1 CoUisionless shock heating of the ions 

Ion heating mechanisms in coUisionless shocks depend on the shock Alfven Mach num- 
ber, the magnetic field inclination angle (&n), plasma parameter (3 and the composition 
of the incoming plasma fiow. The structure of a supercritical shock is governed by the 
ion flows instabilities (see e.g. iKennel et al]|l985l : [Lembege et al.|[2004l : iBurgess et al] 
[2005 '). In a quasi-parallel shock {6n ^ 45°) a mixed effect of a sizeable backstreaming 
ion fraction and the ions scattered by the strong magnetic fleld fluctuations (filling the 
wide shock transition region) results in the heating of ions in the downstream region. 
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The ions reflected and slowed down by an electric potential jump 6(j) at the shock ramp 
of a quasi-perpendicular [On > 45°) shock constitute a multi-stream distribution just 
behind a relatively thin magnetic ramp as it is seen in Fig. [l]and Fig. [2] (left panel). 
The O VII phase densities and distribution functions were simulated with a hybrid code 
for a quasi-perpendicu lar {On = 80°) shock in a hydrogen-helium dominated plasma 
(see iBvkov et aLlbooil - Chapter 7, this volume) . Phase densities x — vx, x~Vy, x^-vz 
of the Ovil ion are shown in Fig. [1] in the reference frame where the particle reflecting 
wall (at far right) is at rest and the shock is moving. The shock is propagating along 
the X -axis from the left to the right and the magnetic field is in the x-z plane. The 
system is periodic in the y dimension. The incoming plasma beam in the simulation 
was composed of protons (90 %), alpha particles (9.9 %) and a dynamically insignifi- 
cant fraction of oxygen ions (O vil) with the upstream plasma parameter /3 ~ 1. The 
ions do not change their initial charge states in a few gyro-periods while crossing the 
cosmological shock ramp where the Coulomb interactions are negligible. 

The simulated data in Fig. [T] show the ion velocities phase mixing resulting in 
a thermal-like broad ion distribution at a distance of some hundreds of ion inertial 
lengths in the shock downstream (see the right panel in Fig. [2]). It is also clear in Fig.[l] 
that the shocked ion distribution tends to have anisotropy of the effective temperature. 
The temperature anisotropy Tj_ ^ 3 Tji relative to the magnetic field was found in that 
simulation. Moreover, the hybrid simulation shows that the T± of the O vil is about 
25 times higher than the effective perpendicular temperature of the protons. Thus the 
ion downstream temperature declines from the linear dependence on the ion mass. The 
si mulations sh o w exc essive heating of heavy ions in comparison with protons. 

iLee fc Wul | 2OO0I I proposed a simplified analytical model to estimate the ion per- 
pendicular temperature dependence on Z/A, where m; — Arup. Specifically, the model 
predicts the ratio of the ion gyration velocity Vig2 in the downstream of a perpendicular 
shock {9n ~ 90°) to the velocity of the incident ion in the shock upstream, vi, 



VI 



A J Bt2 



(3) 



where a = 2eS(j} /mi)V^ < 1, and the potential jump 5(j> is calculated in the shock 
normal frame (see lLee fc WrJ[2000l ). The model is valid for the ions with gyroradii larger 
than the shock transition width Zivi- It is not a fair approximation for the protons, 
but it is much better for temperature estimation of heavy ions just behind the shock 
magnetic ramp. The model of ion heating in the fast, supercritical quasi-perpendicular 
(6*11 > 45°) shocks of Ma > 3 predicts a higher downstream perpendicular temperature 
for the ions with larger A/Z. 



2.2 CoUisionless heating of the electrons 

The initial electron temperature just behind the viscous ion jump of a cosmological 
shock depends on the coUisionless heating of the electrons. The only direct measure- 
ments of the electron heating by coUisionless shocks are those in the Heliosphere. The 
interplanetary shock data compiled by Schwartz ct al. (1988) show a modest, though 
systematic departure of the electron heating from that which would result from the 
approximately constant ratio of the perpendicular temperature to the magnetic field 
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Fig. 1 O VII phase density in a hybrid simulated quasi-perpendicular shock (80° inclination). 
The shock is moving from left to right in the reference frame where the particle reflecting wall 
is at rest. The figures show the oxygen phase densities in x — Vx , x — Vy and x — projections 
from top to bottom respectively. The size of the simulation box in x-dimension is about 300 
li. 



strength (i.e. adiabatic heating). Thus, some modest non-adiabatic electron collision- 
less heating is likely present. In the case of a nonradiative supernova shock propagating 
through partially ionised interstellar medium the ratio Tc/Ti in a thin layer (typically 
< 10^^ cm) just behind a sh ock can be test ed using the optical diagnostics of broad 
and narrow Balmer lines ( e . g . I Ravmond|[200 ll ) . High resolution Hubble Space Telescope 
(HST) Supernova remant (SNR) images ma ke that approach rather attractive. A sim- 
ple scaling Tc/T\ oc was suggested by iGhavamian et aL 1 l|2007l ) to be consistent 
with the optical observations of SNRs. 

Strong shocks are thought to transfer a sizeable fraction of the bulk kinetic energy 
of the flow into large amplitude nonlinear waves in the magnetic ramp region. The ther- 
mal electron velocities in the ambient medium are higher than the shock speed if the 
shock Mach number Ms < \J mp/mc, allowing for a nearly-isotropic angular distribu- 
tion of the electrons. Non-resonant interactions of these electrons with large-amplitude 
turbulent fluctuations in the shock tr ansition region could result in collisio nless heating 
and pre-acceleration of the electrons ( Bvkov fc Uvarovlll999l : lBvkovll2005l ). They calcu- 
lated the electron energy spectrum in the vicinity of the shock waves and showed that 
the heating and pre-acceleration of the electrons occur on a scale of the order of sev- 
eral hundred ion inertial lengths in the vicinity of the viscous discontinuity. Although 
the electron distribution function is in a significantly non-equilibrium state near the 
shock front, its low energy part can be approximated by a Maxwellian distribution. 
The effective electron temperature just behind the front, obtained in this manner, in- 
creases with the shock wave velocity as To oc w^j^ with b < 2. They also showed that if 
the electron transport in the shock transition region is due to turbulent advection by 
strong vortex fluctuations of the scale of about the ion inertial length, then the nonres- 
onant electron heating is rather slow (i.e. b < 0.5). The highly developed vortex-type 
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Fig. 2 Hybrid simulated Ovil distribution function (normalised) as a function of a random 
velocity component 5vy = Vy— < Vy > transverse to the downstream magnetic field in a 
quasi-perpendicular shock (80° inclination). The shock propagates along the x-axis, while the 
initial regular magnetic field is in the x—z plane. In the left panel the distribution in the viscous 
velocity jump is shown. The right panel shows the distribution behind the jump at the position 
of the left end in Fig. [T] Multi-velocity structure of the flow is clearly seen in the left panel, 
while it is relaxing to quasi-Maxwellian in the right panel. 



turbulence is expected to be present in the transition regions of very strong shocks. 
That would imply that the initial Tc/Ti oc v^'^ just behind the transition region 
would decrease with the shock velocity for A^s 3> 1. Here the index a is defined by the 
relation Tj cx v^j^ for a strong shock. The degree of electron-ion equilibration in a coUi- 
sionless shock is a declining function of shock speed. In the case of strong vortex-type 
turbulence in the shock transition region one expects in the standard ion heating case 
with a — 2 and rather small b to have {a — b) < 2. Th at Te/T^ scaling i s some what 
flatter, but roughly consistent with, that advocated by iGhavamian et aL I l|200'if ). On 



the other hand in a coUisionless shock of a moderate strength Ms < W the electron 
transport through the magnetic ramp region could be diffusive, rather than by the 
turbulent advection by a strong vortexes. That results in a larger degree of the col- 
lisionless electron heating /equili bration in the shocks as it is shown in F ig. 4 of the 
paper bv lBvkov fc Uvarov ( 19991 ). Recently. iMarkevitch fc Vikhliniiil ( 2007h argued for 



the coUisionless heating/equilibration of the electron temperature in the bow shock of 
TVfs ~ 3 in the IE 0657 - 56 cluster. 

If the local Mach number A^s of the incoming flow in a strong shock wave exceeds 
y^mp/rric, which could occur in the cluster accretion shocks, the thermal electron dis- 
tribution becomes highly anisotr opic and h i gh fre quency whistler type mode generation 
effects could become important. iLevinsoiil (|l996l ) performed a detailed study of reso- 
nant electron acceleration by the whistler mode for fast MHD shock waves. Electron 
heating and Coulomb rela xation in the strong accretion shocks in clusters of galaxies 



fieatmg and Uoulomb rela xation m the strong 
was discussed in details bv lFox fc Loebl (|l997l ) 
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We summarise this section concluding that a collisionless shock produces in the 
downstream flow a highly non-equilibrium plasma state with strongly different tem- 
peratures of the electrons and ions of different species. Moreover, the postshock ion 
temperatures are anisotropic. The width of the collisionless shock transition region is 
smaller by many orders of magnitude than the Coulomb mean free path (that is of a 
kiloparsec range). We consider now the structure and the processes in the postshock 
Coulomb equilibration layers in the WHIM. 



3 Coulomb relaxation of temperatures in the WHIM 

3.1 Relaxation of the ion temperature anisotropy 

A plasma flow partly randomised just behind the viscous ion jump in a collisionless 
shock transition has an anisotropic velocity distribution with respect to the mean mag- 
netic field. In Fig. [l]the oxygen phase densities in the Vx — x and Vy — x projections 
transverse to the mean field have a substantially wider distribution than that in the 
projection Vz — x parallel to the magnetic field. In many cases the ion velocity dis- 
tributions (like that in Fig. [2] right) can be approximated with a quasi-Maxwellian 
distribution introducing some effective kinetic temperature (more exactly it is the sec- 
ond moment of the distribution). Then the parallel and perpendicular (to the mean 
field) temperatures are different and we approximate the 3D particle distribution as 

/(-X--||)=(2^) (2^) exp(^-^-^y (4) 
We will meas ure the temperatures in energy units in most of the eq uations below 



We will meas ure tne temperatures m energy units m most ot tne eq uations below 
(thus kgT — > T). Ilchimaru fc RosenbluthI 1 1970h (see also a comment bv lKaiserll 19791 ) 



obtained the following equation to describe the ion anisotropy relaxation 

At 2 dt Ti ■ ^ ' 

Here 

_i Syrl/^niZ^e^lnyl 



, t.^l/2„3/2 

15mV T^^ 



(6) 



where In A is the Coulomb logarithm and the effective ion temperature is defined as 



4 ^[(l-^2)y^+^2r||]3/2- 

These equations were obtained under the assumption that the electrons have no dy- 
namical role, but provide a static dielectric background to the ions. We can directly 
apply Eqs. [5l-[7]to isotropisation of the plasma field particles (i.e. protons in our case). 
Isotropisation of the minor ion components is mainly due to their interactions with 
protons and helium because of a low metal number density in cosmic plasma. 
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3.2 Relaxation of the ion and electron temperatures 

Following shock heating and the temperature isotropisation a quasi-Maxwellian dis- 
tribution will be established within all plasma components, i.e. groups of identical 
particles, after a time scale given by Eq. |6l The effective temperatures differs strongly 
between the components. All the plasma particles will undergo Coulomb collisions with 
the protons, alpha particles an d electrons dom inating the WHIM plasma, resulting in 
the temperature equilibration. ISpitzeij (1962) found that the temperature relaxation 



of a test particle of type a with plasma field particles can be approximately estimated 
from 

where 



6 niarrip 

In the thermal equilibrium state the postshock plasma must have a single equilib- 
rium temperature Teq. In a fully ionised plasma without energy exchange with external 
components (i.e. radiation or plasma wave dissipation) Tcq can be found from the con- 
dition of constant pressure in the plane shock downstream resulting in 

Tcq = ^naTao/^?1.a (10) 

In cosmic plasmas it is often a fair app r oxima tion to estimate Tcq from the equation 
2 Tcq — Tc + Tp. Then following ISivukhinI ( 19661 ) the charged particle equilibration can 



be approximately described through the equation 



In 



_t_ 2 / Tn \ „ / To ^ 



+ Cc, (11) 



where Cc is a constant to be determined from the initial temperature Tao of a relaxing 
component a = e,p, 

,^16(2.)V^ npe^ 

Eg. inl and Ea. ll2l allow to calculate the structure of relaxation layers to be seen behind 
a coUisionless shock in the WHIM. In Fig. Owe illustrate the postshock equilibration 
of initially cold electrons with the protons initially heated at the ion viscous jump 
of a coUisionless shock transition. The width of the ion viscous jump in cosmological 
shocks is negligible compared with the equilibration length a^cq = it2''oq, where U2 is the 
downstream flow velocity in the shock rest frame. The characteristic column density 
A'eq = n2Xeci to be traversed by protons and electrons in the downstream plasma (of a 
number density 712) before the temperature equilibration, as it follows from Eg. 1121 does 
not depend on the plasma number density, and A^oq oc w^j^. The corresponding shocked 
WHIM column density can be expressed through the shock velocity vj measured in 
100 kms~^, assuming a strong shock where U2 = 

Acq « 5 X 10" 1-7 [InA]^'^ cm"^, 

or through the shocked WHIM gas temperature Tg (measured in 10^ K): 

Acq ^ 2.5 X 10^^ r|[lnyl]"^ cm~^. 
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Fig. 3 Postshock temperature equilibration between the ion and electron components due to 
the Coulomb interactions. 



The Coulomb logarithm for the WHIM condition is In/l ~ 40. It follows from 
Fig. [3] that in the Coulomb relaxation model the postshock plasma column density 
> SA'oq ensures the equilibration at a level better than 1 %. 

The metal ions can be in itially heated at the shock magnetic ramp to high enough 
temperatures > ATp (see e.g. lKorreck et ahlboOTh for a recent analysis of interplanetary 
coUisionless shock observations with Advanced Composition Explorer). However, in a 
typical case the depth A'^ > SA'^cq is enough for the ion temperature equilibration. 

To study UV and X-ray spectra of the weak systems (Aijjj < 10^^'^ cm^^ and 
-^Ovil ^ ^^^^ cm^^) modelling of shocked filaments of A''ij < 10^^ cm^^ would 
require an account of non-equilibrium effects of low electron temperature Tc/Tcq < 1- 

3.3 Effect of postshock plasma micro-turbulence on the line widths 

We consider above only the WHIM temperature evolution due to the Coulomb equi- 
libration processes. Shocks producing the WHIM could propagate through inhomo- 
geneous (e.g. clumpy) matter. The interaction of a shock with the density inhomo- 
geneities results in the gen eration of MHD-waves (Alfven and magnetosonic) in the 
shock downstream (see e.g. IVainshtein et al.lll99"3 ). The MHD-wave dissipation in the 
shock downstream could selectively heat ions, being a cause of a non-equilibrium To/T; 
ratio. In case of a strong coUisionless shock propagating in a turbulent medium, cos- 
mic ray acc eleration could generate a spect rum of str ong MHD-fluctuations (see e.g. 
iBlandford fc Eichlcr 1987; Bell 2004; Vladi mirov et ak ,2006). These MHD-fiuctuations 
could carry a substantial fraction of the shock ram pressure to the upstream andthen to 
the downstream providing a heating source throughout the downstream. The velocity 
fluctuations could also produce non-thermal broadening of the lines. The amplitude 
of bulk velocity fluctuations is about the Alfven velocity since SB ~ B in the shock 
precursor. The Doppler parameter b derived from high resolution UV spectra of the 
WHIM (see e.g. iLehner et allboOTi : iRichter et al.ll200^ - Chapter 3, this volume): 
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would have in the micro-turbulent limit a non-thermal contribution 

Ont = ~ 



3 6n{p) 5' 

The factor Cv here accounts for the amplitude and spectral shape of the turbulence. For 
a strong MHD-turbulence that was found in th e recent models of stron g coUisionless 
shocks with efficient particle acceleration (e.g. IVladimirov et ahl bood ) one can get 
Ci. ~ 1 in the WHIM, and then 

bnt «4B_9(r2b/i^/0.02)"^/^(l + 2)"^/^(5"^/^ kms"^ 

The estimation of bnt given above for a strong Alfven turbulence may be regarded 
as an upper limit for a system with a modest MHD turbulence. The non-thermal 
Doppler parameter does not depend on the ion mass, but bnt oc S~'^^'^ B{z). Thus 
the account of bnt could be important for high resolution spectroscopy of metal lines, 
especially if the strong shocks can indeed highly amplify local magnetic fields. That 
is still to be confirmed, but a recent high resolution observation of a strong Balmer- 
dominated shock on the eastern side of Tycho's supernova remnant with the Subaru 
Telescope supports the existence of a cosmic ray sh ock precursor where gas is heated 
and accelerated ahead of the shock I Lee et al]|2007l ). High resolution UV spectroscopy 



of the WHIM could allow to constrain the intergalactic magnetic field. Internal shocks 
in hot X-ray clusters of galaxies have modest Mach numbers and the effect of the 
Alfven turbulence is likely less prominent than that in the strong accretion shocks in 
clusters and in the cosmic web filaments. X-ray line broadening by large sca l e bulk 
motions in the hot intracluster medium was discussed in detail bv lFox fc Loebl (|l997l ) 
and llnogamov fc Sunvaev, (|2003i ') . 



4 lonisation state of the WHIM 

To s imulate absorption spe ctra of bright quasars in the intervening WHIM filaments 
(e.g. iKawahara et al.l (|2006h as an example of such a modelling) one should solve the 
ionisation balance equations for the charge states of metal ions with ac count taken of 
all the LTE processes and also the nonthermal particle contribution (see lPorquet et al] 
( 200lh for a discussion of a role of non-relativistic super-thermal distributions). The 
ionisation balance equation can be written as 

riq = rlc [nq_lCq_l — WqCq — UqOq + nq+lOq+l] -|- 

\ ^ r i/ion /Trrcc , Trion\ , t ^rcc 1 , 

+ '^j |_'^q-l^j,q~-l -"q(^j.q +^,q )+'^q+l^j,q+lJ + 

j=H,Hc,Ho+ 
+nq_li?q_l — UqRq. 

Here q is the charge state of an element, Cq is the coUisional (and autoionisation) rate 
q q + I, Qq are the radiative and the dielectronic ionisation rates g — > g — 1 (in 
cm'' s^^), Vj'^q'^. '^j'q" are the charge exchange rates with the ion j (in cm'^ s~^) and -Rq 
is the photoionisation rate of an ion (in s~^). T he rates of different processes can be 
calculated for different temperature regimes (see iKaastra et al]|2008l - Chapter 9, this 
volume). We just limit our discussion here to one example of such a simulation. 

In Fig.|4]we illustrate coUisional ionisation equilibrium curves of oxygen ions in the 
present epoch (at z = 0) as a function of WHIM density S. The various boundaries 
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1 + 6 

Fig. 4 Phase diagram for oxygon in the IGM at redshift z = 0. Density is parameterised by 
the S value as was defined above. The soUd hne in the upper left hand corner labelled 'fd < tn' 
indicates where electron and proton fiuids reach equilibrium in a Hubble time (tn). Low density 
gas will not radiatively cool over a Hubble time to the left of the boundaries marked 'tcool > tn'- 
The three curves are labelled with the metallicity, Z, expressed as a fraction of Solar metallicity. 
The solid curve labelled ' adiabatic evolution' indicates the locus of gas that has only undergone 
adiabatic compression or expansion since high redshift (initial condition Te ^ 10^ K). Critical 
boundaries for the ionisation equilibrium of oxygen are: the shaded area in the upper right 
hand corner indicates the regime where the coUisional ionisation timescale is shorter than the 
photoionisation timescale, for ionisation O VIII — > Oix. The two boundaries labelled 'tcoii < ^h' 
indicate where the coUisional ionisation timescale becomes shorter than the Hubble time. At 
lower temperature, the ionisation balance cannot be in (coUisional) equilibrium. Upper (solid) 
curve is for ionisation Ovil (0+^) — ► Ovill (0+^), lower (dashed) curve for Ovi (0+^) — ► 
Ovil (0+^). The steep solid curve labelled 'troc (Oxygen) > tjj' indicates where the radiative 
recombination timescale (O IX — > O Vlll) exceeds the Hubble time (no recombination at low 
densities). 



separate regimes under which a certain process does, or does not attain equilibrium 
over a Hubble time. We show in Fig.|4]some critical boundaries for kinetic and thermal 
equilibrium. The solid line in the upper left hand corner labelled 'toi < tjj' indicates 
where electron and proton fluids reach kinetic equilibrium (proton temperature equal 
to electron temperature) in a Hubble time (tjj): at low density and high temperature, 
such equilibrium does not attain. Low density gas will not radiatively cool over a 
Hubble time to the left of the boundaries marked 'tcool > ^h'- The cooling time was 
calculated for coUisionally ionised gas. The solid curve labelled ^adiabatic evolution' 
indicates the locus of gas that has only undergone adiabatic compression or expansion 
since high redshift (initial condition To ~ K); all shock-heated gas will be above 
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this lino right after passing through a shock. The shaded area in the upper right hand 
corner indicates the regime where the colhsional ionisation timescale is shorter than the 
photoionisation timescale, for ionisation Oviii — > Oix. The two boundaries labelled 
'^coU < *h' indicate where the colhsional ionisation timescale becomes shorter than the 
Hubble time. At lower temperature, the ionisation balance cannot be in (colhsional) 
equilibrium. . The steep solid curve labelled 'tree (Oxygen) > tjj' indicates where the 
radiative recombination timescale (O ix — > O Viii) exceeds the Hubble time since there 
is no recombination at low densities. 

Note that wc illustrate here only the colhsional equilibrium case. For the more 
appropriate case of radiative cooling in photoionisation equilibrium, the cooling times 
will be even longer, due to the fact that a photoionised pleisma is highly overionised 
compared to the characteristic ionisation- and excitation potentials, which suppresses 
the (very effective) colhsional cooling contribution. In our figure, the lines tcooj = tn will 
shift to the right if we calculate with the probably more realistic case of photoionisation 
equilibrium. However, as we argued above the colhsional equilibrium is also of interest, 
since it represents a conservative case. The range of ionisation states of oxygen in the 
WHIM filaments can be observed in the absorption spectra of bright quasars. 

5 Conclusions 

We discussed some specific features of the WHIM heating processes by collisionless 
plasma shocks driven by gravitationally accelerated flows in the LSS structure forma- 
tion scenario. 

• Collisionless shock heating of ions in the WHIM results in a highly non-equilibrium 
initial state with a strongly anisotropic quasi-Maxwcllian ion distribution just behind 
a very thin magnetic ramp region. The ion temperatures just behind the shock depend 
on the ion atomic weight, the charge state and the shock magnetic field inclination. 
The ion temperatures could decline from a simple linear scaling with the ion mass 
providing an excessive heating of heavy ions. 

• The ion and electron temperatures relax to the equilibrium state through Coulomb 
collisions in a layer of the depth A^oq ~ lO^^v^ cm~^ behind a shock of velocity wy. 

• Strong collisionless plasma shocks with an efficient Fermi acceleration of energetic 
particles can generate strong MHD waves in the downstream region that will result in 
a non-thermal broadening of the emission/absorption lines in the observed WHIM 
spectra. 
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